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We speculate about the possible interpretations of the recently ob¬ 
served excess in the h ^ Tfi decay. We derive a robust lower bound on 
the Higgs boson coupling strength to a tan and a muon, even in presence 
of the most general new physics affecting other Higgs properties. Then 
we reevaluate complementary indirect constraints coming from low energy 
observables as well as from theoretical considerations. In particular, the 
tentative signal should lead to r —)■ /ly at rates which could be observed 
at Belle H. In turn we show that, barring fine-tuned cancellations, the 
effect can be accommodated within models with an extended scalar sec¬ 
tor. These general conclusions are demonstrated in explicit new physics 
models. Finally we show how, given the h —)■ r/i signal, the current and 
future searches for /i —)■ ey and /x —>■ e nuclear conversions unambiguously 
constrain the allowed rates for h —)■ re. 
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1 Introduction 


The discovery of the Higgs boson 0i offers numerous observables where the validity 
of the Standard Model (SM) can be tested. The CMS collaboration has recently 
reported a slight excess with a significance of 2.4 cr in the search for LFV decay 
h —)■ r/i [^. The best fit for the branching ratio of the Higgs boson to rp under 
assumption of SM Higgs production is found to be 

e(ft^T,,) = (o.84is|?)%. (1) 

This recent hint has received great attention in the literature (see and references 
therein). An ATLAS study of h —)■ r/r in the hadronic r channel observes no 
excess and is consistent with the CMS result ([^. 

In light of the tentative signal it is instructive to revisit compatibility of such large 
i3(h —)• r/i) with other low-energy lepton flavor violation (LFV) probes in the context 
of popular New physics (NP) models. 

2 Higgs effective theory approach 

The mass terms and Higgs boson couplings of charged leptons after electroweak sym¬ 
metry breaking (EWSB) can be parametrized in the general case as 

Cf = -mA/A - Vij (AC) h + ...+ h.c., (2) 

where the ellipsis denotes non-renormalizable terns involving more than one Higgs 
boson. While in the SM the Yukawas ?/ ~ m are diagonal in the mass basis, NP 
could misalign the two matrices and via and/or induce h ^ ry decays with a 
branching ratio of B{h —)■ r/i) = (IW/iP + \yA)- assumes that the total 

Higgs boson decay width (F/i) is SM-like and enlarged only by the contribution from 
h —r/i then the measurement in Eq. (|T|) can be interpreted as a two-sided bound, 
0.0019(0.0008) < AAFAhA < 0 . 0032 ( 0 . 0035 ) at 68% (95%) C.L.. In general, 
the experimentally measured h ^ ry branching fraction depends also on other Higgs 
couplings contributing both to its total decay width (Th) as well as its production 
cross-section {(Th)- In particular, a given signal can be reproduced for larger (smaller) 
values of \yr^l\ and \yfMr\ by enhancing (suppressing) Th and/or suppressing (enhancing) 
ah- Individual effects of F^ and ah can be disentangled by performing a global fit to 
all Higgs production and decay event yields at the LHC |^. Numerically, we find in 
this case 

0.0017(0.0007) < < 0.0036(0.0047) at 68% (95%) C.L.. (3) 

In the following we assume the SM contains all the relevant degrees of freedom 
at energies (P(few 100) GeV, whereas for the additional heavy degrees of freedom we 
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assume they been integrated out. The natural ranges for the effective Higgs couplings 
follow from the hierarchy between the muon and tan lepton masses , 

^J\yr^.y,r\ ^ =0.0018. (4) 

It is interesting to note that almost the whole parameter space in allowed by h —)■ rp 
is also compatible with the above naturalness criterium [^. 

Phenomenologically the most relevant low-energy constraints on this scenario are 
the one- and two-loop contributions to the operators LH{a-B)E and LTaH{a■W°‘)E, 
where and are the hypercharge and weak isospin held 

strengths, respectively, and are the Pauli matrices. These operators can mediate 
the strongly constrained radiative LFV decays. The most stringent constraint comes 
from the r —)• /ly decay [^, mediated by the effective Lagrangian 

= ClQl'^i + CrQr-i + h..C. ^ (5) 

where Ql, R^ = {e/87r‘^)mr{fi(j'^^PL,Rr)Eap, Pl,r = (1 T75)/2 and Eajs is the electro¬ 
magnetic held strength tensor. The resulting EFT correlation between B{h — rfi) in 
and B{t —)■ py) is shown in left-hand panel in Fig. (diagonal dashed orange line), 
assuming SM values of all Higgs boson couplings except and y^r- In fhe same 
plot, the CMS preferred range of B{h —)■ rp) in Eq. ([^ is displayed by the horizontal 
blue band, while the current (B(t —)■ py) < 4.4 x 10“® @ 90% C.L.) and projected 
future {B{t —)■ py) < 3 x 10“® @ 90% C.L.) indirect constraints are shaded in 
light and dark gray vertical bands, respectively. We observe that within the EFT 
approach, the CMS signal is well compatible with the non-observation of r —)■ py at 
the B factories and will marginally remain so even at Belle H. 


3 Type-Ill Two Higgs Doublet Model 


Extensions of the SM with an additional SU{2)l doublet scalar are an ehective de¬ 
scription of several NP scenarios (e.g. supersymmetric extensions or models address¬ 
ing the strong CP problem via a Peccei-Quinn symmetry), for a recent review c.f. (^. 
The most extensively studied version of the Two Higgs Doublet Model (THDM) is 
the type-H model in which one of the doublets couples to up-type quarks, while the 
other one couples to down-type quarks and charged leptons, avoiding in this way the 
tree level flavor changing neutral currents (FCNCs). All LFV in this case comes at 
1-loop and is aligned with the SM flavor structure thus rendering LFV processes to 
a negligible level, compared to the one of the SM [^. In type-HI THDM Yukawa 
couplings are generic and allow for large LFV. In the case with a MSSM-like Higgs 


potential the following relations hold 11 
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Figure 1: Left-hand side panel: Correlation between B{h — )■ rp) and B{t — )■ /X7) in 
various NP scenarios. The present experimental result for B{h ^ Tfi) is shown in 
horizontal blue band [^. Current and future projections for B{t ^ n'j) experimental 
sensitivity are represented with vertical light and dark gray bands, respectively. 
Superimposed are the predictions within the EFT approach (diagonal dashed orange 
line), in the type-III THDM (green and black bands), and in models with scalar 
leptoquarks (diagonal red and orange shaded band). Right-hand side panel: Allowed 
region in the B{h — )■ Te)-B{h —)■ Tfi) plane when experimental upper bounds on 
/i —)■ ey and /i — e conversion rates are taken into account. Pink region is permitted 
in the effective theory setting while the dashed line indicates how much the region 
will shrink if Mu2e and MEG II experiments see no signal events. Green region is 
allowed within type-III THDM model with niA = 0.5 TeV and tan/3 = 10. Rulers 
indicate how much the region shrinks with increasing tan f3 or while dashed lines 
correspond to improved experimental upper bounds on /i —)■ ey and /i — e as described 
in the text. 


while for the masses 


^ +m‘^w 


and m% = m\ -\- m| — m 


— iit^ -r iit^ — Here /3 is the 
angle that diagonalizes the mass matrices of charged scalars and pseudoscalars while 
a is an analogous angle for the neutral scalars. The relevant part for the discussion 
of LFV are the Yukawa couplings in the charged lepton mass eigenstate basis 
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where the LFV Yukawa couplings can be written as 
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The off-diagonal parameters drive the LFV phenomena, while the coefficients 
for Hk = {H, h, A) can be found in |^. Using the above relations, we find the tree-level 
result for the h —)■ r/i branching fraction, 

rrih 


B{h —>• r/i) = 


IbTrTft 


(sin a tan/3-f cos (|( 
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We have found that even for light pseudoscalar masses {mA), modihcations of hWW 
and hZZ relative to their SM values is negligible |^. 

The decay width of r —>■ py is driven by a one- and two-loop Barr-Zee diagrams 
with virtual charged or neutral Higgses. At one-loop the amplitude is suppressed by 
two small Yukawa couplings, while the two-loop result is proportional to yu and a 
single LFV Yukawa. Indeed one hnds that Barr-Zee contributions dominate the rate 
of r —)■ j^. We sample the parameter space of ef.,. for a chosen values 

of tan/3 and m^. The ranges allowed for required to fulhll the naturalness 

criterium of Eq. (|^ and stay within the perturbative regime. 

The 1-loop amplitude of the r —)■ /ly process further depends on the diagonal 
Yukawa coupling which is experimentally constrained by the searches for h ^ tt 
decays by the ATLAS 


12 and CMS 13 


experiments. The naively averaged signal 
strength of two rr signal strengths results in = 1 . 021 q 2 o which directly constrains 
y^^. A scenario with SM-like y^^ coupling corresponding to = 1 for hxed tan/3 = 
10 and niA = 0.3 TeV is presented in left-hand side panel in Fig. by a black narrow 
stripe. This scenario easily passes both experimental constraints. On the other hand, 
for masses signihcantly larger than 500 GeV it is not possible to reconcile both 
predictions with the corresponding experimental values. Allowing to departure 
from zero within the allowed range we obtain correlation represented by the green 
band in left-hand side panel in Fig. [T] In particular, this additional freedom in r —)■ py 
breaks the strict correlation with the h —)■ rp rate. 


4 Scalar Leptoquarks 

A scalar leptoquark state (LQ) can induce h —>■ r/i decay via quark-LQ penguin 
diagrams. Inspection of the helicity structure of the relevant diagrams reveals that 
both chiralities of leptons and top quark have to couple to the LQ state, in order keep 
the leptoquark couplings perturbative 0. LQ state can couple to the Higgs also via 
the “Higgs portal” operator, that comes with an additional parameter, 

A. 

4.1 The Ai = (3,1, —1/3) case 

The Yukawa couplings of Ai are given by the following Lagrangian 

U, = + i/gr AiU^^' + h.c., (10) 
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where Q* = (Wi, and ?7* = are the quark weak doublets and up-quark singlets, 
respectively. We explicitly show flavor indices i, j = 1,2,3, and SU(2) indices a,b = 
1, 2, with ei 2 = 1. The h ^ Tfi decay width in the presence of Ai scalar is then, 

T{h rp) = {\y^y^r\^ + btrl/jr) • (11) 


Here the relevant loop function further depends on the portal coupling A and is given 
in [^. The state Ai also contributes to the r —)■ pq via same y couplings as the ones 
present in h —)■ rp: 


B{t fi-f) 
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where Xt = ml/m\_^ and expression for the hi function can be found in Ref. j^. 

The impact of a non-zero Higgs portal coupling A in the scenario with Ai has also 
been studied. As an example, for tuai = 650 GeV the loop function dependence on 
the portal coupling is gi = —(0.26 -|- 0.12A). Thus, a positive large A could relax the 
leptoquark Yukawa couplings and yield sizable h —)• r/i rates without violating the 
r —>■ /iq constraint. However, the Higgs portal coupling also induces corrections to 
the h —)■ qq decay and to gluon-gluon fusion (ggF) induced Higgs production with the 
leptoquark running in the triangular loop. Taking those modifications of the Higgs 
couplings into account we have found that one can compensate small leptoquark 
Yukawas by scaling up A, which is itself bounded from above only by perturbativity 
requirements j^. 

The correlation between the h —)■ r/i and r —>■ /iq branching ratios in presence 
of the (3,1,—1/3) leptoquark state for |A| < 1 and ttiai > 600 GeV is depicted in 
left-hand side panel in Fig. [^with a pink stripe, demonstrating that r —)■ /iq basically 
excludes this LQ state as an explanation of h —)■ r/x signal. 


4.2 The A 2 = (3, 2,7/6) case 

The Yukawa couplings of the A 2 leptoquark to SM fermions are 

Aa, = - yjr A“6“'L^’' + h.c., (13) 

The h —)■ rp decay rate in this case is 

F(h ^ rp) = {\y^^tyfr? + b/tl/jr) • (14) 

Allowed values for the Higgs portal coupling A can be inferred from a global fit to the 
Higgs data as has been done for the portal coupling of the (3,1, —1/3) state. 
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The decay width of r —)■ /ly are proportional to the couplings responsible for 
h — r/i: 


^ h7) 


am^ 





( 15 ) 


In this leptoquark scenario the bound on —)■ /ly) excludes sizable B{h —)■ rp) 

due to the strict correlation between the two observables. See the orange stripe in 
left-hand side panel in Fig. where the portal coupling is restricted to |A| < 1. 


5 h ^ Tfi VS. h ^ re 


A positive indication for h ^ rji decay can be combined with stringent experimental 
limits on /i—e LFV processes to constrain r—e processes. In particular, in models with 
tree-level h —)■ r/x (EFT, THDM III) the product of the B{h —)■ re) and B{h —)■ r/i) 
is bounded from above by the rates of /i —)■ ey and fi — e conversion on nuclei. This is 
due to the fact that tree level Higgs decays to r/i (re) depend on (l/er,Te) while 

the same sets of couplings contribute at the loop level to /i —)■ ey and /x — e conversion 
via diagrams with a virtual r. In the effective theory framework, the contributions 
to the /i —)■ ey process stemming from a virtual r are m,- enhanced with respect to 
diagrams with intermediate /x or e states leading to 


C rs-* 

r — 


nir —3 -|- iXr — — 2 log Xj 


m. 


8m^(l — XrY 


V/XTyr 
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where we have neglected the effects of the light lepton masses. The coefficient is 
obtained from Eq. (16) by replacing yij —)■ y*^. The /x —>■ ey branching fraction is thus 
sensitive to a distinct combination of the LFV Yukawas: 


H(/X ^ ey) ~ {ly.rVrel^ + \yr,yer\^) , = 185 . (17) 


On the other hand, /x — e conversion on nuclei is most sensitive to vector cur¬ 
rent effective operators (ey^PL^^/x) The branching fraction i3(/i —e)Au = 

F(/i -)■ e)Au/rcapture Au cau be put in the form 


H(/i ^ e)Au = B^^ (|//er2/p.P + \yreyr^x\") , B^^ = 4.67 X lO"" , (18) 


where the relevant numerics have been taken from Ref. (also c.f. |^). The com¬ 
plementary information on the LFV couplings extracted from /x —)■ ey and /x — e 
conversion allows for correlation between the Higgs LFV h —)■ r/x and h —)■ re decays: 


B{h —)■ Ty)B{h —)■ re) 


8 X 10“^° 


H(/i ^ ey) 
10-13 


+ 3 X 10“^ 


H(/i —)■ e)Au 

KFi3 


(19) 
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The best experimental limit on —)■ e)Au < 7 x 10“^^ (at 90% C.L.) was achieved 

by the SINDRUM II Collaboration while the best upper bound on i3(/i —)■ ey) < 
5.7 X 10“^^ (at 90% C.L.) was recently determined by the MEG Collaboration (l^ . 
Note that with the current experimental data the sum on the right-hand side of 
Eq. (19) is completely saturated by the fi — e conversion contribution. Combining 
the two bounds with the central value for the h ^ TjjL branching fraction in Eq. ([^ 
leads to an upper bound B{h —)■ re) < 0.26. This is above the current indirect 
constraint coming from searches for r —>■ eq which reads B{h —)■ re) < 0.19 . Future 
improvements of bounds on /i —>■ eq and especially fi — e, and assuming B{h —)■ Tfi) 
stays at the percent level, will indirectly probe B{h —re) at the 10“® level [^. Similar 
analysis can be carried over to the THDM III framework, where the resulting bound 
on h —re renders this decay invisible at the LHC (green region in right-hand side 
panel in Fig. [^. 


6 Conclusions 

Motivated by the experimental indication of h —)■ r/i events by the CMS Collaboration 
we have examined the implications of LEV Higgs decays at the percent level on 
several extensions of the SM. We have shown how a tentative B{h —)■ Tfi) signal can 
be combined with other Higgs measurements to yield a robust lower bound on the 
effective LEV Higgs Yukawa couplings to taus and muons. In explicit NP models, 
the r —)■ pq constraint is generically more severe. In fact, an eventual observation of 
h —)■ r/i at the LHC together with indirect constraints would point in the direction 
of an extended SM scalar sector, minimal example being THDM of type HI. We 
have also examined models where h —)■ r/i is generated at loop level to demonstrate 
difficulties with these models. Finally, we have combined the signal of h —)■ rp with 
experimental limits on p —)■ eq decays and fi — e conversions in nuclei to yield robust 
bounds on B{h —re). In particular, considering only the Higgs EFT effects, the two 
LEV Higgs decay rates could still be comparable. On the other hand, the THDM HI 
cannot accommodate both branching ratios at the percent level. It turns out that 
currently planned improvements in experimental searches for fi — e LEV processes 
will be able to probe the product B{h —)■ Tfi)B{h —)■ re) at the 10“^ level in generic 
EFT and to order 10“^^ or better within the THDM HI. 
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